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We present an extension of the continuous damage fiber bundle model to describe the gradual 
degradation of highly heterogeneous materials under an increasing external load. Breaking of a fiber 
in the model is preceded by a sequence of partial failure events occurring at random threshold values. 
In order to capture the subsequent propagation and arrest of cracks, furthermore, the disorder of 
the number of degradation steps of material constituents, the failure thresholds of single fibers are 
sorted into ascending order and their total number is a Poissonian distributed random variable over 
the fibers. Analytical and numerical calculations showed that the failure process of the system is 
governed by extreme value statistics, which has a substantial effect on the macroscopic constitutive 
behaviour and on the microscopic bursting activity as well. 

PACS numbers: 46.50+a, 62.20.Mk, 81.40.Np 
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I. INTRODUCTION 

Materials with a highly disordered microstructure ex- 
hibit a variety of constitutive characteristics when sub- 
jected to an increasing external load — from perfectly 
brittle to perfectly plastic as well as strain hardening and 
softening [U, Experiments have revealed that even if 
the constituents are brittle, the macroscopic behaviour 
shows these variations, especially if the internal struc- 
ture of the specimen has a hierarchy of length scales. 
A possible explanation of this observation is the grad- 
ual activation of failure mechanisms relevant at different 
length scales of the materials' microstructure [Sj. This 
interaction of microstructure and failure mechanism is 
typical for fiber reinforced composites, where fibers of a 
few micrometer diameter are embedded in a matrix ma- 
terial Q • The same geometrical arrangement is found on 
a larger scale when pieces of timber are glued together to 
form high strength wood structures [5j. In such materi- 
als the microscopic origin of gradual degradation is the 
accumulation of damage in fibers and timberpieces due 
to the growing population of microcracks [6|, |7| , further- 
more, the growth and arrest of larger cracks spanning 
over several length scales of the specimen . 

Fiber bundle models (FBM) are one of the most ade- 
quate approaches to understand the fracture of heteroge- 
nous materials [u!, [H, [H, [H. Recently a continuous 
damage fiber bundle model (CDFBM) has been intro- 
duced where the stiffness of fibers is reduced in subse- 
quent failure events representin g th e gradual degradation 
mechanism of the material [3. llSl. [l6j . Recent experi- 
ments revealed that due to the highly heterogeneous mi- 
crostructure of faults, defects (like the knotholes of wood, 
[13]), and also particularly failure resistant spots, the 



maximum number of failure of composites' constituents 
show strong variations inside samples, which has also 
an important effect on the evolution and propagation of 
cracks. Since the stress concentration around the tip of 
the crack increases with the crack length, cracks that 
grow larger than a certain critical size are unstable 
However, cracks can also become arrested at strong lo- 
cations inside the material, where the stress concentra- 
tion around the growing crack demands an increasing 
strength of subsequent arrest locations [13, [3 ■ This 
phenomenon therefore requires the knowledge of the dis- 
tribution of the strongest points in the sample, an issue 
that is related to order statistics [19, uO| . 
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In this paper we propose an extension of CDFBM in or- 
der to capture the experimentally relevant failure mech- 
anisms of composite materials discussed above. We con- 
sider a bundle of parallel fibers which undergo a sequence 
of failure events when the local load on them exceeds 
the corresponding threshold values. The maximum num- 
ber of allowed failures is a Poisson distributed random 
variable in the model which accounts for the disorder of 
degradation steps of fibers in composites. Since grow- 
ing cracks sweep over locations of increasing strength, 
in the model the breaking thresholds of consecutive fail- 
ure events of single fibers are sorted in ascending order. 
We explore the complexity of the model by analytical 
and numerical calculations and show that the extreme 
value statistics of failure thresholds has a substantial ef- 
fect on the fracture process. Contrary to the conven- 
tional CDFBM, the macroscopic response of the system 
does not have a plastic plateau, but instead strain hard- 
ening occurs. Computer simulations revealed that on 
the microlevel the burst size distribution follows a power 
law with a crossover of the exponent from 5/2 to 3/2 
when changing the parameters of the model. Beyond the 
general interest, our model can in particular be applied 
to describe the stability of large scale composite wood 
structures. Following the catastrophic collapse of sev- 
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eral buildings composed of such structures in the winter 
2006/2007, see e. g. [2l|, this issue has received growing 



interest in the recent past. 



II. RANDOMLY DISTRIBUTED FAILURE 
NUMBERS 



As a model we employ a continuous damage fiber bun- 
dle model (CDFBM), in which the damage law of the 
classical dry fiber bundle model (DFBM) is supplemented 
by a gradual degradation of fiber strength in the sequence 
of failure events [1, [l^ . It was shown in Ref . [T^ that for 
certain choices of the model parameters a variety of ex- 
perimental situations can be recovered, i.e. either strain 
hardening or plasticity can occur. On the micro-scale, 
the size distribution of avalanche events shows a power 
law behaviour, but the exponent is different from the or- 
dinary FBM [24I, and for certain choices of parameters 
an exponential cutoff appears . 

The CDFBM is constructed as follows: the bundle con- 
sists of N parallel fibers on a square lattice with identi- 
cal Young- modulus E and random failure thresholds crj^, 
i = 1^ . . . , N with a probability density p and distribu- 
tion function P. Under loading, the fibers behave linearly 
elastic until they reach their respective points of failure 
and break in a brittle manner, i.e. as soon as the load on 
a fiber exceeds its breaking threshold cr^^, the fiber will 
fail. If the external strain is kept fixed (strain-controlled 
loading) , then the fibers break one by one in the order of 
their breaking thresholds, and the full constitutive curve 
is explored. Under stress controlled loading, however, the 
load on a breaking fiber is redistributed on the remaining 
intact ones, where it can trigger subsequent avalanches of 
fiber breakings [23*| . At the maximum of the constitutive 
curve, a catastrophic avalanche then results in failure of 
the entire system [l^, [23| . 

The failure law of the DFBM is modified in the 
CDFBM by assuming that at the failure point the Young 
modulus of the fiber is reduced by a factor a, where 
< a < 1; consequently the stiffness of the fiber af- 
ter failure is aE. The loading of the fiber will then re- 
sume in a linear manner with the reduced stiffness until 
the next breaking threshold is reached. The parameter 
^max determines the maximum number of failures allowed 
for a single fiber. The damage threshold crj^ can either 
be kept constant for all the breakings (quenched disor- 
der) or new failure thresholds of the same distribution 
can be chosen (annealed disorder) after each instant of 
failure, which can model a microscopic rearrangement of 
the material after failure, cf. cases (a) and (b) of Fig. [T] 

It can be assumed that in an actual experimental sit- 
uation, the number of times that a constituent of the 
material can break is an independent realization of an 
integer random variable. A prime example is the frac- 
ture of wood, specifically of glued timber d, |33| , where 
only a few large defects and the finite number of glued 
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FIG. 1: Damage law for a single fiber for the cases of quenched 
(a), annealed (b) disorder, cf. 0], and for the model with 
sorted thresholds (c) discussed in Sec. IIIII In all cases, af- 
ter breaking of a fiber upon reaching a failure threshold ctj^, 
loading is resumed with a stifi'ness that is reduced by a factor 
a. 



joints determine the extreme value statistics that gov- 
erns the propagation and arrest of cracks. This fact can 
be incorporated by modeling /cmax as a random number, 
which is governed by a Poissonian distribution 



(1) 



A new parameter then enters the model, which is the 
mean value of fcmax, = (fcmax)- 

With this prescription, the macroscopic response of the 
system can be expressed as 



a = e 



CXD u _^ 



fe=0 



(2) 



i. e. the Poissonian distribution is convoluted with the 
degradation term obtained for the constitutive behaviour 
in the continuous damage model @. In the case of 
quenched disorder for the failure thresholds, the prob- 
abilities Pki^) that at a given deformation s a randomly 
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chosen fiber has failed exactly k times is 



I -Pie) 
Pk{e) = { P(a'=-i£) - P(a%) 



,fc = 0; 
,l<k<k„ 
h — k 



(3) 



For the following discussion, the damage thresholds will 
be drawn exclusively from a WeibuU distribution with 
A = l,m = 2, unless otherwise mentioned. 

We can also apply this model to the case of annealed 
disorder, where 



[1 - Pia'^e)] nto ^(a'e) , < A: < fc,nax - 1; 



n 



J=0 



P{a^e) 



, /c — k„ 



(4) 

however, we restrict this discussion to the quenched dis- 
order case. In Sec. IIIIl a new concept will be introduced 
which to a certain extent describes a third alternative to 
the cases of quenched and annealed disorder. 

In the analytical solution for the constitutive be- 
haviour, Eq. (12), two physically strongly distinct cases 
can be realized by appropriate choices of the summation 
limit of the innermost, bracketed term: if the summa- 
tion extends from zero to fcmaxj as indicated, the fibers 
will retain a residual stiffness after the limiting case of 
^max failures, i.e. hardening of the fiber bundle occurs in 
the limit of large e, and the asymptotic behaviour of the 
bundle is described by 



^asympt. 
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= £e-"(i-"^ (5) 

Fig. [5] demonstrates the perfect agreement between the 
analytical solution, Eq. ([2]), and a strain controlled 
numerical simulation, where the asymptotic linear be- 
haviour of Eq. ([5]) is recovered. In order to model the 
failure of materials, however, the failure law has to be 
slightly modified: after k* = fc,„ax — 1 failures, the load 
on a fiber must be set to zero, and the constitutive be- 
haviour changes to 



a = £ 



.fc=0 



(6) 



Again, in this case the constitutive curves displayed in 
Fig. [5] show an excellent agreement between the analyt- 
ical solution Eq. ([6]) and the simulation data. At first 
glance, the hardening behaviour that emerges after the 
main course of loading may appear to contradict the fact 
that residual stiffness is not explicitly taken into account 
and the terms Pfc(e) with k — /cmax are excluded in the 
failure law, Eq. ^\ this regime is dominated by the fibers 
with fcinax = 0, i. e. , fibers that never break, and since the 
expectation value of these fibers 

the asymptotic behaviour even in the case without an 
explicit residual stiffness reads 
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FIG. 2: Constitutive behaviour for a fiber bundle with 
quenched disorder and Poisson distributed fcmax for three val- 
ues of K, with (r.s.) and without (n.r.s.) residual stifi^ness. 
Symbols: analytical solution, Eq. ((2)|, solid lines: simulation 
results. The dashed lines show the asymptotic behaviour in 
the presence of residual stifi'ness, Eq. ([Sjl, and without, Eq.((7|. 



Apparently, the dominance of the fibers with fcmax = 
diminishes with increasing k, and at least for this choice 
of the disorder distribution, under stress controlled load- 
ing the hardening regime cannot be accessed, i.e. all 
fibers break before traversing the local minima of the 
slope. One should note that — with or without a residual 
stiffness term — the fibers with vanishing fc^ax can be 
excluded from both the simulations and the analytical 
calculations, in case of which the hardening behaviour 
in the second case will also disappear. The distribution 
«K (fcmax) will not be purely Poissonian anymore, and an 
account of this situation is planned for a future publica- 
tion [m. 

Concerning the avalanche size distribution, this model 
reproduces the behaviour of the case of fixed fcmax, which 
has been discussed in There, it was found that 

for larger values of a, corresponding to a > 0.3 for 
the WeibuU distribution, the distribution can be fitted 
to a power law, the exponent of which also depends on 



For small values of fcr, 



the usual mean field be- 



asympt 



ee 



(7) 



haviour with an exponent —5/2 is obtained [22| , whereas 
for larger values of fcmax a smaller exponent w 2.12 ap- 
pears. Examining the effect of the Poissonian term with 
a choice of k = (fcmax) corresponding to the ordinary 
CDFBM, we find in Fig. [3] a quantitative agreement be- 
tween the two models, i. e. the Poissonian term causes 
no visible change to the avalanche statistics, and even 
the crossover between the two power law exponents is 
recovered. 
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FIG. 3: Avalanche size distribution for a fiber bundle of size 
L — 201 with Poisson distributed fcmax for different values 
of K. Inset: Avalanche size distribution for a fiber bundle of 
size L = 201 with fixed fcmax, where the fcmax values shown 
correspond to the k values of the main plot. 



FIG. 4: Marginal densities for the order statistics of a WeibuU 
distribution with \ — l,m = 2, for the respective values 
j = 1, . . . ,n — 4 (solid lines), together with the underlying 
PDF (dashed line), cf. Eq. @. 



III. SORTED FAILURE THRESHOLDS 

In highly disordered materials subjected to an increas- 
ing external load, cracks nucleate in the early stages of 
loading at the weakest locations in a spatially random 
manner. As the load increases, simultaneously to the 
nucleation of new microcracks the existing cracks prop- 
agate and become unstable. Advancing cracks can be 
arrested by high strength locations of the material. Be- 
fore macroscopic failure occurs, advancing cracks can un- 
dergo several activation and arrest events. Since stress 
concentration at the crack tip increases as the crack be- 
comes longer jjy|, arresting can only be realized by local 
material elements of increasing strength. These growth 
and arrest events result in a gradual degradation of the 
macroscopic sample strength. 

In order to provide a more realistic representation 
of this gradual degradation process sweeping through 
material elements in the increasing order of their local 
strength, we further extend the CDFBM by sorting the 
activation thresholds into increasing order. Fig. [1] case 
(c), provides a graphical illustration of the sorting and 
the ensuing damage law for a single fiber. It is important 
to emphasize that from a physical point of view this case 
is a mixture of the annealed and quenched disorder cases 
discussed previously 0]. On the one hand, the sorted 
model bears resemblance to annealed disorder, since the 
consecutive thresholds are different from each other; on 
the other hand, it could also be classified as quenched as 
the thresholds are fixed in advance. Sorting of a series 
of random numbers imposes a correlation between these 
numbers, and we will have to resort to a mathematical 



theorem from the field of order statistics in order to ob- 
tain a complete understanding of the failure mechanism. 
It should be stressed that in the following discussion we 
will consider fcmax to assume a fixed value, although the 
addition of the randomly distributed fcmax can be triv- 
ially incorporated. Here a single fiber — seen as a meta- 
element representing smaller constituents — models the 
progress of a crack, which does not proceed continuously 
but comes to a halt at certain values of the fibers' strain 
e. We may apply the damage law of the CDFBM, but im- 
pose the additional condition that the load on the fiber at 
subsequent instants of the arrest should increase. Hence 
we can draw the failure thresholds from a random distri- 
bution, for which we will use again the WeibuU distribu- 
tion, and store them in sorted order. As in the previous 
discussion, the WeibuU distribution employed will have 
the parameters A = 1 and to = 2. 



A. Macroscopic Response 

As mentioned above, bringing an array of n random 
numbers in sorted order necessarily invokes correlations 
between them, and the distribution of the random variant 
at the ith position is not governed by the PDF of the 
unsorted random numbers anymore. 

From order statistics, the following results pertaining 
to a sequence of random numbers Xi , . . . , Xn drawn from 
a CDF F{x) and density f{x) are known [2^: if the 
Xi are brought into increasing order, then the numbers 
< ■ • ■ < X(^n) are called the order statistics of the 
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FIG. 5: Constitutive curves for the model with sorted random 
thresholds, without residual stiffness. Solid lines: simulation 
results, symbols: analytical solution Eq. (|12p . dashed and 
dotted lines: simulations for the conventional CDFBM with- 
out sorting, shown for comparison. The failure thresholds are 
drawn from a WeibuU distribution with A = 1, m = 2. 



Xi, which have the joint density 



1=1 



gix(i), ■ ■ • ,a;(„)) = < 



0, 



'CO < X^i) < ■ ■ ■ < .X(„) < CO 



elsewhere . 



(8) 

The marginal density for the jth order statistics ^(j), 
1 < j < n is 



— OO < t < CO . 



[F{t)y-'[i~F{t)r-\f{t), 



(9) 



The marginal statistics Eq. @ is therefore the ade- 
quate replacement of the PDF, i. e. the distribution func- 
tion for the random number at the j-th position, if n 
random numbers have been drawn. In order to illus- 
trate this result, in Fig. |4] the marginal statistics for the 
jth random number, 1 < j < n = A, together with the 
underlying WeibuU distribution with A = 1, m = 2 is 
shown. It is apparent from Fig. |4] that with increasing 
j, the marginal statistics share the peaked characteris- 
tics of the underlying PDF, and that the position of the 
maxima reflects the sorting. Also, the marginal statistics 
become wider with increasing j, although this effect is 
not too pronounced. 

Having found an analytical expression for the marginal 
statistics, the constitutive behaviour of this model can be 
expressed in a closed form. In analogy to the annealed 



case of the CDFBM, we denote by Pfc(e) the probability 
that a fiber has failed exactly k times at a strain e: 



Pk{e) = I 



[l-G(fe+i)(a^-e)] ri'G(,+i)(a%), 

< fc < fciiiax ^ 1 

fcmax-1 

n G(j+i)(a^e), 

h — h 



(10) 



where 



Gq-)(x) 



gj (t) dt 



(11) 



is the integral associated with the marginal statistics 
9{j)W^ corresponding to the CDF of unordered random 
numbers. The second case in Eq. (jlOp with k = fcmax cor- 
responds to the residual stiffness of the bundle. With this 
result the analytic solution for the constitutive behaviour 
reads 



ais) 



fcmax -1 

E 

A:=0 



G(fc+i)(a'' 



fe-i 

J=0 



(12) 

if the hardening term is skipped in order to account for 
material failure. 

It should be noted, though, that in the formula for 
the fc-th failure probability Pfe(e), Eq. PU)) . the integral 
Eq. pT|) appears, which cannot in general be solved an- 
alytically due to the structure of the integrand, Eq. 
In order to obtain the constitutive behaviour the inte- 
gral Eq. (fTTjl has to be evaluated numerically. In Fig. [5l 
the stress-strain curves obtained in this way are plotted 
for three values of fcmax, where an excellent agreement 
between the numerical and analytical results, Eq. 
can be observed. It can be seen in Fi g, p that con- 
trary to the case of unsorted thresholds [l5[ — whether 
quenched or annealed disorder is of no importance — that 
the constitutive curve does not develop a plateau, it al- 
ways increases monotonically and has a quadratic maxi- 
mum where macroscopic failure occurs. As a consequence 
of extreme value statistics, with growing fcmax the criti- 
cal stress iTc and strain Ec increase, indicating a higher 
macroscopic load bearing capacity. It can be seen from 
the general expression of the constitutive curve, Eq. (|12p . 
that the macroscopic failure of the system is mainly con- 
trolled by the largest thresholds whose distribution can 
be obtained from Eq. ([9]), setting j = fcmax- Analyzing 
the constitutive behaviour of the system considering only 
the largest thresholds, j = fcmax, yields 



aX [ln(fcn 



1)] 



1/m 



and 



(13) 



(14) 



for the failure stress and strain, respectively, assuming 
WeibuU distributed failure thresholds with parameters 
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FIG. 6: The critical stress ctc and strain Ec for a bundle with 
sorted random thresholds from a WeibuU distribution with 
A = l,m = 2 and without residual stiffness as a function of 
the maximum number fcmax of failures. The symbols represent 
simulation results, whereas the lines display evaluations of 
Eqs. (|13I14[) and numerically obtained maxima of Eq. (|15|) 
with the setting k = fcmax, respectively. 
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FIG. 7: Constitutive curve (solid line) and avalanche sizes 
without residual stiffness as a function of strain e from a sim- 
ulation with fcmax = 60, m — 2. Inset: ratio Ei+i/si for sub- 
sequent bursts with Ai > 200, for the values a = 0.7 (unfilled 
circles) and a — 0.8 (filled squares), where the lines indicate 
the respective values of 1/a. 



A, m. Eq. ^ implies that the sorted CDFBM does 
not have a plastic limit as in the conventional CDFBM, 
i. e. no plateau of the a{e) emerges. Instead, the strength 
of the bundle is an asymptotically increasing function of 
fcmax, namely, CTc increases logarithmically whereas £c in- 
creases exponentially with fcmax- This is illustrated in 



b 




FIG. 8: (Color online) Constitutive curve for the model with 
sorted failure thresholds and residual stiffness with A = 1, 
m = 2 for various values of fcmax in a double-logarithmic 
plot, under stress controlled loading. The curves display the 
same asymptotic behaviour as the residual stiffness curves in 
Fig. [21 although it is not visible in this logarithmic represen- 
tation. The inset shows a section of the constitutive curves 
on a logarithmic-linear scale, where the onset of oscillations 
with increasing fcmax is better visible. 



Fig. [6l where the values Uc and Sc obtained by com- 
puter simulations are compared to the analytical results, 
Eqs. (fT3TT4l) . 

For very high values of fcmax, a distinguished regime of 
oscillations appears in the constitutive curves, see Fig. [T] 
This plot shows the constitutive behaviour in a stress 
controlled simulation with fcmax — 60, together with the 
avalanche sizes A recorded at each loading step. Appar- 
ently, the constitutive curve displays a large amount of 
oscillations with horizontal plateaus, which coincide with 
large scale bursts of breaking events. The position of the 
peaks suggests a regularity of some kind. In order to 
quantify this regularity of the peak events, the inset car- 
ries information about the ratios Ei+i/si of subsequent 
failure events A with A^ > 2000. Obviously, this ratio 
assumes a constant value of « 1/a after a brief onset pe- 
riod, where a = 0.8 is the load reduction parameter used 
in all the simulations presented in this discussion, and a 
comparison with the case of a = 0.7 is presented in order 
to confirm the influence of the load reduction parameter. 
It should also be noted that the envelope of the constitu- 
tive curve remains monotonically increasing. In order to 
analyze this oscillation phenomenon an investigation on 
the evolution of the breakdown process is necessary. 

In fact, an analytical argument can be made about 
the origin of these oscillations, and their position relative 
to the state of loading e. In the constitutive formula, 
Eq. (|12p . there appears a product of integral marginal 
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distributions Gq+i) (a-'e), where < j < fc — 1, and 
< fc < /Cinax — 1- From the structure of Eq. and 
if a WeibuU distribution with m = 2, A = 1 is assumed 
as usual, it can be understood that the respective deriva- 
tives (x) take on their maximum vahie even for very 
large j, n, and j n at a numerical value of the argument 
X = 0(1). Therefore, the g{j+i) for large values of j pos- 
sess well defined peaks at very large values of e, where 
a% w 1. The positions of these peaks become strongly 
separated for subsequent indices j, j -I- 1 if j = 0{n) and 
n is large. Consequently, in Eq. the product of the 
integral quantities 0(^+1) (a-'e) can be replaced by the 
largest factor Gfj.) (a'^^^e) for relatively large k, and to- 
gether with the leading term e[l — G(fc+i) (a'^e)] a peak 
structure is formed. It was confirmed numerically for 
fcmax = 60 and k > 30 that the positions e of these max- 
ima, which are defined through the function 

rukie) = a''e[l - G^k+i){^''e)]G^k){a''~'e) (15) 

practically coincide with the observed peaks of 
avalanches, furthermore, setting k = fcmax also yields 
good estimates for the critical stress and strain, see Fig.|6l 
An oscillatory regime also appears if residual stiffness 
is present, as shown in Fig. [5] for stress controlled simu- 
lations for various values of fcmax, where oscillations are 
clearly visible for fcmax ^ 30, as in the case without resid- 
ual stiffness. Again, the appearance of oscillations de- 
pends on the choice of fcmax, and sets in at fcmax ~ 30. 
With increasing fcmax also pronounced strain hardening 
occurs, whereas in the conventional CDFBM a plastic 
plateau is present. After passing the strain hardening 
regime, the bundles attain an asymptotic regime of con- 
stant slope a'^'"""' for all values of fcmax: where also macro- 
scopic failure occurs for all cases investigated. It should 
be noted that therefore at low values of fcmax, i- e. without 
oscillations, the slope of the constitutive curve is always 
finite positive, whereas for higher values of fcmax, there 
are sections of the constitutive curve with zero slope due 
to the oscillations. 



B. Bursts of Fiber Breakings 

The presence of the oscillations, and therefore of a lo- 
cally vanishing slope, has a distinguished effect on the 
distribution of avalanche sizes, as demonstrated by Fig. [9] 
for the case with residual stiffness, and in Fig. [TO] for 
the case without; this is also suggested by the presence 
of large size avalanches in Fig. [7] (no residual stiffness). 
It can be seen that in the presence of residual stiffness, 
Fig. El two remarkable features are present. First, for 
all values of fcmax, there is a regime with a power law of 
exponent —3/2 for small avalanche sizes. Secondly, for 
small values of fcmax, an exponential cutoff appears for the 
larger avalanches; the presence of the initial power law 
regime with a following cutoff is confirmed by rescaling 
both axes by the average size of the largest avalanche. 




FIG. 9: Avalanche size distributions for the CDFBM with 
sorted failure thresholds for various values of fcmax, 1 < 
fcmax < 60 and residual stiffness. Results from stress con- 
trolled simulations of L — 201 fiber bundles, averaged over 
100 realizations. 



(Amax). It is important to note that the distributions 
D{A) can be collapsed onto a single master curve, which 
was then fit by the functional form 

D{A) cx (j^) " exp L^—) , (16) 

V (Amax) / \ C (Amax) / 

as can be seen in Fig. Illi where the curves with low val- 
ues fcmax < 20 have been used. The above arguments are 
supported by the fact that the power law of exponent 
3/2 with an exponential cutoff provides a perfect fit to 
the master curve obtained numerically. It has to be em- 
phasized that in the fitting, the value of the exponent was 
fixed to 3/2, and the fit was obtained solely by varying 
the parameter c in Eq. (jl6p . 

For high values of fcmax, i. e. in the presence of oscilla- 
tions, a crossover is observed from an initial regime with 
a power law of exponent —3/2, to another regime with 
the mean field power law of exponent —5/2, and finally 
to a peaked regime for very large avalanches of about the 
system size. 

A similarly complex behaviour is observed in the cases 
without residual stiffness, see Fig. [TO] There, for low 
values of fcmax, the usual mean field behaviour of a power 
law with an exponent —5/2 is observed [13, [24|. For 
intermediate values of fcmax, a crossover occurs between 
an initial —3/2 power law part to a mean field part for 
larger avalanches, and the position of the crossover shifts 
to larger avalanche sizes with increasing fcmax. For higher 
values of fcmax, a limiting curve with a crossover at A « 
10^ can be identified, and again a peak of avalanches of 
the order of the system size is found for fcmax = 60. 
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FIG. 10: Avalanche size distributions for the CDFBM with 
sorted failure thresholds and no residual stiffness for various 
values of fcmax, 1 < fcmax < 60. Results from stress controlled 
simulations of L = 201 fiber bundles, averaged over 200 real- 
izations. 
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FIG. 11: Rescaled avalanche size distributions for the 
CDFBM with sorted failure thresholds and residual stiffness, 
for small values of fcmax outside the oscillatory regime. A best 
fit with the numerical values a = 1.22 and l3 = 1.5 for the 
exponents has been used to obtain the collapse of the curves 
and the dashed fit curve. 



The features found in the avalanche size distributions 
can all be explained on the basis of the fine structure of 
the oscillations, which is displayed in Fig. [T^ for both 
stress and strain controlled loading, with fcmax = 15 and 
the choice m — 10 for the WeibuU parameter. From 
this illustration it becomes apparent that the horizon- 
tal plateaus in the stress controlled simulations actually 
correspond to regions of decreasing stress a under strain 
controlled loading, which cannot be accessed in the stress 
controlled mode. Also, for these parameter values, the 
integral Eq. (fTTjl could be solved analytically by using a 
computer algebra program, and the analytic solution of 
Eq. (fT2|l shown in Fig. [12] is in excellent agreement with 
the simulation results. 

In general, a fiber bundle model can only produce large 
avalanches if the constitutive curve has at least one max- 
imum, where the susceptibility to a small increment of 
the external force diverges [12, . Avalanches with a 
power law distribution are generated in the vicinity of 
the maximum of cr(e), where the shape of the maximum 
determines the value of the exponent r. Quadratic max- 
ima typically result in r — 5/2, the value obtained in 
the absence of both oscillations and a residual stiffness 
term. If loading is stopped at a strain before reaching 
the maximum, i. e. before global failure occurs, an expo- 
nential cutoff in the avalanche size distribution appears 
[T^ l. which is visible in Figs. l^landfTTl where due to the 
residual stiffness term the bundle fails macroscopically 
after passing exclusively through regions of finite posi- 
tive slope, without quadratic maxima, as for the cases 
fcmax < 30 no oscillations occur. 

In the oscillatory regime, however, a{e) passes a se- 



ries of consecutive maxima with an increasing amplitude. 
Under stress controlled loading the system jumps from a 
local maximum of (j{e) to the ascending side of the next 
maximum which is somewhat higher than the previous 
one, see Fig. 1121 The jump implies that a large amount 
of fiber breakings occur in a single avalanche removing 
all fibers which have breaking thresholds lower than the 
load of the ending point of the jump on the next peak 
of cr(e). Consequently, when loading is continued along 
the ascending side of the peak determined by Eq. (|15p . 
the response of the system is determined by a disorder 
distribution which is critical in the sense of Refs. [13, HI] , 
i. e. weak fibers are removed so that the lower cutoff of 
the disorder distribution falls close to the local critical 
deformation, the location of the next peak. As it has 
been shown in (s^ . |35| . when the disorder distribution 
approaches criticality, the avalanche size distribution ex- 
hibits a crossover from a power law with an exponent 
r = 3/2 for the small avalanches, to another exponent of 
T — 5/2 for the large ones, irrespective of the effective 
range of interaction. This effect can be recognized in the 
fcmax = 60 curve in Fig. [Q] 

The same argumentation holds for the case without 
residual stiffness: here, a global quadratic maximum is al- 
ways present, so for low values of fcmax the mean field ex- 
ponent —5/2 is found. However, in the presence of oscil- 
latory structures, there appear a series of local quadratic 
maxima, and the stepping effect described above yields 
a series of local critical threshold distributions, which re- 
sults in the crossover of exponents visible in the avalanche 
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FIG. 12: Section of the constitutive curve for the model with 
fcmax = 15, m — 10 and no residual stiffness. Solid line: 
strain controlled simulation, dashed line: stress controlled 
simulation, dots: exact evaluation of the analytical solution, 
Eq. (Ull). 



size distribution, cf. Fig. [TOl 

It has to be stressed that the effects found here in the 
CDFBM with sorting have no equivalent counterparts in 
the conventional CDFBM without sorting. There, the 
macroscopic behaviour yields a plastic plateau and no 
steps appear in the constitutive curve; consequently, the 
size distribution of avalanches shows no signatures of crit- 
icality [lB|. The existence of oscillations, synchronized 
avalanche bursts and a critical crossover of the avalanche 
size distribution exponents is a genuine peculiarity of the 
order statistics accompanying sorting. 

The analytical argumentation presented above also 
makes it clear that oscillations cannot appear for all dis- 
order distributions and all values of fcmax; rather, its 
appearance is restricted to combinations of either large 
A^max and high disorder corresponding to low values of the 
WeibuU parameter m, or smaller values of fcmax and cor- 
respondingly low disorder, such that the maxima defined 
through Eq. ([TS]) are clearly separated. Fig, [represents 
a numerical survey of the apparent occurrence of oscil- 
lations. One can see that a well defined and smooth 
separatrix can be found which isolates the regime with 
oscillations from the regime without in the {/cmax, "ti} pa- 
rameter space. 

IV. CONCLUSIONS 

Motivated by experimental observations on the frac- 
ture process of composite systems having a hierarchy of 
length scales, we extended the continuous damage fiber 
bundle model by taking into account that a hierarchy of 



FIG. 13: Phase diagram denoting the appearance of oscilla- 
tions in the constitutive curves of the model with sorting and 
no residual stiffness, and the presence of pronounced spikes 
in the avalanche size distribution. Filled circles: oscillatory 
regime; crosses: normal regime; unfilled circles: regimes can- 
not be distinguished. The solid line denotes the approximate 
location of the separatrix between the two regimes. 



length and energy scales for damage initiation and crack 
growth and propagation exists. Since macroscopic dam- 
age and arrest event are based on large defects and high 
strength zones in the material, their number is not a con- 
stant but varies between samples of the same production 
batch. Therefore, two new features have been added to 
the classical CDFBM, and their effect on the microscopic 
and the macroscopic damage evolution has been investi- 
gated. 

First, the maximum number of failures fcmax has been 
modeled as a Poissonian random variable, which incor- 
porates the existence of disorder with respect to the fi- 
nite number of macroscopic defects and strength zones 
which can effectively govern the macroscopic breaking of 
certain materials. The presence of the Poissonian term 
has a distinct effect on the constitutive behaviour as it 
induces a hardening regime, which becomes more pro- 
nounced for small average values of fcmax- For the mi- 
croscopic behaviour, the introduction of the Poissonian 
distributed fcmax leaves the distribution of avalanche sizes 
invariant, and a crossover from a power law with an expo- 
nent —5/2 to a power law with another exponent —2.12 
for increasing values of < fcmax > is found, in analogy to 
the conventional CDFBM. 

In a next step we introduced sorted failure thresholds 
in the model and explored the inherently complex be- 
haviour in this case by analytical and numerical means. 
A parameter regime has been identified where the dam- 
age evolution of all fibers synchronizes and considerable 
changes to the microscopic quantities can be observed, 
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depending on the amount of disorder and the maximum 
number of aUowed failures. It was shown that the ex- 
treme value statistics of failure thresholds has a substan- 
tial effect on the fracture process of the system both on 
the micro- and macro-level. This theoretical study can be 
particularly important for composite materials produced 
by assembling components with a large variation in their 
respective defects. The application of this model to the 
failure process of large wooden structures is in progress. 
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